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We use first-principles density functional theory (DFT) calculations to investigate the ground 
state structures of both BaCeOa (BC) and Pd-doped BC (BCP). The relaxed structures match 
closely with recent experimental scattering studies, and also provide a local picture of how the 
BC perovskite lattice accommodates Pd. Both stoichiometric and oxygen- deficient materials are 
considered, and structures with an O vacancy adjacent to each Pd are predicted to be favored. 
The oxidation state of Pd in each doped structure is investigated through a structural analysis, 
the results of which are supported by an orbital-resolved projected density of states. The vacancy 
stabilization by Pd in BCP is explained through redox chemistry and lattice strain relief. 



I. INTRODUCTION 



The structure and composition of solid oxide per- 
ovskites have been varied widely, in order to tailor their 
properties for applications [HEHSHHSIISIEIEIE]. The 
development of perovskites as either catalysts themselves 
or as better oxide supports for catalytic applications is 
the focus of great current interest. [TOl HTJ |T2] 

Nishihata et al. reported [T3l HI] that when precious 
metals (Pt, Rh, or Pd) are mixed into solid oxides of 
La(Fe,Co)03, then sintered, the precious metal becomes 
incorporated within the host perovskite as an oxidized 
species. When subjected to a reducing atmosphere, the 
precious metal is extruded from the host matrix and de- 
posited on the surface of the oxide as a fully reduced 
metal. According to Nishihata et al. this metal is the 
active species for catalytic reactions, such as the oxida- 
tion of CO to CO2. When the reaction is complete, the 
catalyst is subject to an oxidizing atmosphere. This ox- 
idizes the precious metal, which returns back into the 
host oxide, where it can be stored for future catalysis. 
The two-stage cycle can be repeated many times. This in- 
creases the efficiency of their catalyst system, presumably 
because incorporation inhibits precious metal agglomera- 
tion under oxidizing conditions. As a result, the catalytic 
metal retains higher surface area for much longer than on 
conventional catalytic substrates. This work in so-called 
"intelligent catalysts" has focused mainly on La-based 
perovskite hosts, where precious metal dopant levels of 
up to 5% were achieved. 

Based on crystal chemical size rules, we investigated 
whether the upper limit of 5% precious metal doping is 
related to the size mismatch between Pd and Fe. For 
analysis of ionic size in perovskites, the tolerance factor 
t is the central parameter: 



where Ra-o is the sum of A and O ionic radii and Rb-o 
is the sum of B and O ionic radii. LaFeOs (rLa=l-36 A, 
rpe ^0.6 A, ro=1.42 A) has a tolerance factor of 0.97. 
The introduction of a Pd+^ on the B site (rpd=0.86 A) 
would locally make t = 0.86 which would create signifi- 
cant strain and limit solubility. Pd+^ is a smaller ion but 
perhaps a more difficult oxidation state to stabilize. 

To improve on the La-based materials and to find new 
oxide functionality, we explored the ability of other ma- 
terials to incorporate precious metals reversibly. One im- 
portant modification scheme has been the use of large A- 
and B-site cations, to increase the unit cell volume. It 
has been shown that large cations can yield perovskite 
oxides that function well as either proton or oxygen ion 
conductors [I5l|l6], the most promising being doped BC. 
The presence of large ions (rBa=l-61 A and rce=0.87 A) 
on both A- and B-sites makes the host perovskite cell 
unusually large, and we propose that this large volume 
can more readily accommodate precious metal ions. 

Another important perovskite design strategy pursued 
in this work is the use of an ionic, low- valence A-site 
cation (Ba), in order to stabilize higher redox states in 
adjacent cations, in this case the cation incorporation of 
Pd. The large size of BC and the redox-stabilization of 
Ba make BC a suitable host for precious metal dopants. 
We recently investigated [iTl [18] how these design prin- 
ciples apply to Pd-doped BC, namely whether BCP re- 
versibly incorporates Pd and whether it is an active cat- 
alyst. 

BaCeOa has t=0.94 and is predicted to be an anti- 
ferrodistortive dielectric material. Since Ce+^ and Pd+^ 
have the same ionic radii in 6-fold coordination, Pd dop- 
ing maintains t and should not induce much strain. To 



2 



maintain charge neutrality, Pd+^ would need to be ac- 
companied by an oxygen vacancy. 

In this study, we examine the atomic and electronic 
properties of Pd-doped BC in detail using DFT. The local 
bonding environments around Pd as well as the electron 
fillings of atomic orbitals on Pd and its neighbors are ex- 
amined. A particular focus is the character of the highest 
occupied molecular orbital (HOMO) and the few lowest 
unoccupied molecular orbitals (LUMO). These electronic 
states are compared with expectations based on crystal 
field splittings in the computed atomic geometries. This 
combined structural and electronic computational study 
helps explain the Pd incorporation previously reported. 

There are several theoretical studies of pure BaCeOs 
and with various dopants [19, 20j |21| |22J, but none of 
these reported theoretical predictions of the full ground 
state structure and symmetry group of pure BaCeOs us- 
ing DFT. Our preliminary theoretical studies are the first 
to consider Pd doping of BaCeOs, which agree well with 
recently published experimental data [171 IlH]- In this 
paper we investigate the atomic and electronic structure 
of BaCeOa and BaCei-a^Pda^Oa in more depth, building 
upon previously reported results. 



II. 



METHODOLOGY 



In this study, an in-house solid state DFT code em- 
ployed in previous studies [23l |24] is used to relax the 
ionic positions and optimize lattice constants. The gen- 
eralized gradient approximation (GGA) [25 of the ex- 
change correlation functional and a 4 x 4 x 4 Monkhorst- 
Pack sampling of the Brillouin zone [26 are used for 
all calculations. All atoms are represented by norm- 
conserving optimized [27 designed nonlocal [28 pseu- 
dopotentials. All pseudopotentials are generated by the 
OPIUM code [29j. The calculations are performed with 
a plane wave cutoff of 50 Ry. Undoped and Pd-doped 
BaCeOa calculations were performed on forty atom unit 
cells in a 2 X 2 X 2 supercell arrangement allowing for the 
three-dimensional octahedral tilt that is seen in experi- 
ment. In the Pd-doped perovskite, two supercells of com- 
position BaCeo.875Pdo.125O2.875 were studied, as well as 
the structure containing no vacancy, BaCeo.875Pdo.i25 03. 
Structure 1 contains an O vacancy adjacent to Pd, and 
Structure 2 contains the vacancy in the next unit over, 
between two Ce. Structure 3 contains no vacancy. 



III. 



RESULTS 



In the forty-atom BaCeOa 2x2x2 unit cell, the Oe 
tilt angles are all 12.2°. However, the tilts along (100) 
and (001) are anti-phase and indistinguishable from each 
other, while the tilts along (010) are in-phase. This struc- 
ture belongs to the a~a^a~ tilt system. The converged 
structure belongs to the Pnma symmetry group, and the 
forty-atom computational cell we used is a supercell of 





Experiment 


Theory 


a 


6.21 


6.28 


b 


8.77 


8.81 


c 


6.23 


6.30 


Ba(x) 


0.0207 


0.0273 


Ba(z) 


-0.0075 


-0.0077 


01(x) 


-0.0132 


-0.020 


oi(^) 


0.4250 


0.4194 


02(x) 


0.2771 


0.2833 


02(^) 


0.0383 


0.0446 


02{z) 


0.7239 


0.7184 



TABLE I: Comparison of experimental [17] and theoretical 
structural parameters for BaCeOa. 



the primitive ^2 x 2 x ^2 cell. Relaxed DFT lattice 
constants agree well with experiment [17] , as seen in Ta- 
ble [T[ as do the positions of ions from the experimental 
refinements. 

Structures 1, 2, and 3 were modeled with a = 8.799, 
b = 8.759 and c = 8.823 A , the relaxed DFT lattice 
constants of BaCeo.875Pdo.i25 03 (Structure 3). Doping 
BaCeOs with cationic Pd reduces the symmetry of the 
unit cell from Pnma to the space group PI. The tilt 
system in BaCeOa is severely interrupted with the ad- 
dition of one Pd, then becomes closer to the pure BC 
tilt pattern with the addition of a Pd-0 vacancy pair. In 
Structure 3, the substitution of one Ce for one Pd with no 
accompanying vacancy has tilt angles that range from 2 
to 15°. The average tilt angle of this system is 9.6 ± 4.6°. 
This perturbation results in a dramatic decrease in sym- 
metry, also displayed by the Ce-0 bond lengths. The 



three nearest Ce neighbors (Ce4, 6 and 7 in Table IV) all 
have two distinct sets of bond lengths; four bonds around 
2.21 A and two long bonds around 2.36 A . The second 
nearest neighbors have six similar bond lengths of 2.28 A 
on average. The three nearest Ce neighbors to Pd are en- 
closed within octahedra that distort, then tilt, as shown 
in Figure [3| The distorted Ce octahedra are the result 
of the formation of strong, short bonds (~ 2.03 A) to 
the small Pd^+ cation. The bond lengths of 2.03 A are 
typical of Pd^+ in an Oe environment (ro2-=1.42 and 
rpd4+ =0.615). 

Of the two structures with an oxygen vacancy, the 
higher energy structure is Structure 2, where the O va- 
cancy is in between two Ce. This is 0.95 eV higher 
in energy than Structure 1, because the reduction of 
Pd^+ to Pd^+ is more favorable than the reduction of 
two Ce^+. The vacancy in Structure 1 is further stabi- 
lized by the formation of a square planar environment 
around Pd^+. As shown in Figure [l] the Pd in Struc- 
ture 1 is four-coordinate square-planar, and not five- 
coordinate square-pyramidal. In a square-planar envi- 
ronment, rpD2+ = 0.64. This is supported by the Pd-0 
bond lengths in Structure 1 from Table [ll| The first four 
shortest bond lengths are on average 2.04 A, and the 
fifth is 2.66 A. This longer bond length occurs as a re- 
sult of the apical O in the CeOs square pyramid moving 



Cel 
Ce2 
Ce3 
Ce4 
Ce5 
Ce6 
Ce7 



Pdl 



2.26 
2.23 
2.23 
2.22 
2.24 
2.05 
2.22 



2.27 2.27 2.27 

2.24 2.24 2.27 

2.24 2.24 2.25 

2.22 2.22 2.25 

2.24 2.25 2.25 

2.21 2.21 2.23 

2.22 2.22 2.25 



2.27 2.32 
2.27 2.29 

2.26 2.28 
2.35 2.36 

2.27 2.28 
2.23 4.51 
2.35 2.35 



2.04 2.04 2.04 2.05 2.66 4.30 



TABLE II: Six shortest B-0 distances (in A), for Structure 1, 
BaCeo.875Pdo.125O2.875, with vacancy between Pdl and Ce6. 




FIG. 1: Structure 1, where Pd dopant is adjacent to an O 
vacancy. This shows the lesser degree of octahedral distortion 
of the nearest neighbors when compared to Structure 2, as 
well as the increasingly symmetric tilt. Pd is green, Ba are 
light blue, Ce are gray and O red. 



away from the Pd. This does not occur in Structure 2, 
where the vacancy is in between two Ce (Figure [2|. All 
six Pd-0 bond lengths in Structure 2 are around 2.02 A, 
whereas the two CeOs surrounding the vacancy seem to 
contract inward, as seen from the bond lengths in Ta- 



ble III These bond lengths are similar to those seen in 
Structure 3, where there is no vacancy. 

This clearly indicates that Pd reduction via an O va- 



Cel 
Ce2 
Ce3 
Ce4 
Ce5 
Ce6 
Ce7 



Pdl 



2.21 2.22 

2.22 2.22 
2.17 2.21 
2.26 2.28 
2.22 2.22 
2.26 2.27 
2.22 2.22 



2.22 2.22 
2.28 2.28 

2.21 2.22 

2.28 2.31 
2.26 2.27 

2.29 2.31 

2.22 2.22 



2.24 4.34 

2.30 2.34 

2.23 4.34 

2.37 2.38 

2.32 2.35 

2.35 2.37 

2.36 2.42 



2.01 2.01 2.01 2.02 2.02 2.03 



TABLE III: Six shortest B-0 distances (in A), for Structure 
2 BaCeo.875Pdo.125O2.875, with vacancy in between Cel and 
Ce3. 




FIG. 2: Structure 2 has an O vacancy between two Ce atoms, 
not adjacent to Pd (lower plane, not shown). Ba are light 
blue, Ce are gray and O red. 



Cel 
Ce2 
Ce3 
Ce4 
Ce5 
Ce6 
Ce7 



Pdl 



2.28 
2.23 
2.25 
2.20 
2.24 
2.21 
2.21 



2.28 2.28 

2.24 2.26 

2.25 2.25 

2.20 2.22 
2.24 2.26 

2.21 2.21 
2.21 2.22 



2.28 2.29 2.29 

2.26 2.27 2.27 

2.25 2.26 2.26 
2.22 2.36 2.36 

2.26 2.26 2.26 

2.21 2.37 2.37 

2.22 2.35 2.35 



2.03 2.03 2.03 2.03 2.04 2.04 



TABLE IV: Six shortest B-0 distances (in A), in Structure 3 
BaCeo.875Pdo.i2503, with no vacancies 



cancy may be needed to alleviate some of the strain of 
this system induced by the small Pd^+ ion. When an O 
vacancy is placed between two Ce, as in Structure 2, the 
closest corner sharing octahedra decrease their tilt angles 
to either 9 or 2°. The remaining octahedra slightly in- 
crease their tilt angles to either 13 or 14°. Structure 2 
is more distorted than Structure 3, and it seems to have 
less of a concerted tilt preference. The average tilt an- 
gle in Structure 2 is 9.9 ± 5.1°. In Structure 1, when 
the vacancy is between Pd and Ce, the corner sharing 
octahedra are not as affected, decreasing their tilt angles 
by at most 2°. The CeOs next to the vacancy decreases 
its tilt angle to 4°. The average tilt angle of Structure 
1 is 9.9 ± 3.4°. The average tilt angle of Structures 1 
and 2 is the same, but the standard deviation is larger 
for Structure 2. The largest interruption of concerted 
tilt systems by a vacancy occurs in Structure 2, where 
it is not stabilized by the Pd^+ cation. All four super- 
cell structures resemble the tilt system a~a+a~, however 
symmetry has been broken by both the Pd substitution 
and accompanying O vacancy. 

To elucidate the local structure of a disordered per- 
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FIG. 4: Experimental neutron PDF G{r) of BCP, shown 
in dashed blue line, and DFT PDF G{r) of Structure 1, 
shown in solid red line. Experimental BCP was of formula 
BaCei-ccPdcc03-5 where x=0.10. 



FIG. 3: Structure 3, where Pd dopant is not accompanied by 
an O vacancy. Pd is green, Ba are light blue, Ce are gray and 
O red. 



ovskite, such as BCP, we simulate a neutron-scattering 
pair distribution function (PDF) by combining DFT- 
computed atomic positions and known neutron scattering 
factors. To make direct comparison with experimental 
neutron PDF, we subtract the average number density 
and multiply by Anr: 



G{r) = 47rr(p(r) - po{r)) 



(2) 



where p{r) is the microscopic pair density and po{r) is the 
average number density. The experimental PDF G(r) of 
BCP generated with the program PDFFIT [30 matches 
closely the DFT-derived PDF of Structure 1 generated by 
an in- house code used in previous studies [STJ |32] . The 
correspondence in peaks, as shown in Figure [4j is strong, 
especially for short r. This further supports the agree- 
ment between experimental data and our calculations, as 
well as the assignment of Pd^+ being stabilized by a O 
vacancy in Structure 1. The match between experiment 
and theory was much less close for PDFs generated from 
Structure 2 and Structure 3, providing evidence that Pd- 
vacancy pairs are present in the experimentally prepared 
BCP sample. 

Analysis of the density of states (DOS) of the three 
structures shows that Structure 1 contains an interesting 
feature absent from either Structures 2 or 3. Above the 
Fermi level, there is a (i-orbital local ized on Pd that is 
d^2_y2 .This is shown in Figure 5(c) 



and indicates the 



presence of Pd^+. This (i-orbital is not present in either 



the HOMO (Figure |5(a)|) or LUMO (Figure |5(b)|). In a 
square-planar environment, this is the non-degenerate d- 
orbital that is highest in energy. It is absent from both 



Structures 2 and 3 because Pd instead assumes a Pd^+ 
oxidation state in octahedral coordination. 

A comparison of the projected density of states of the 
Pd (i-orbitals in Figure [6] shows that the Pd in Structure 
1 is different than the more similar Structures 2 and 3. 
Since Structures 2 and 3 both contain Pd^+ in an octa- 
hedral environment, they have the same d-states, though 
shifted in energy. The Pd in Structure 2 is slightly more 
ionic than in Structure 3, an effect caused by the O va- 
cancy in Structure 2. This O vacancy, in between two Ce, 
has caused the octahedron around Pd to contract slightly, 
forming stronger, more ionic, Pd-0 bonds. In Structure 
1, new d-states closer to the Fermi level are populated, 
and previously unfilled bands, above the Fermi level, de- 
crease in intensity. This is because Pd has been stabilized 
in a square planar environment as Pd^+. The shift of the 
d-states closer towards the Fermi level when compared to 
Structures 2 and 3 also support Pd being less ionic. 



IV. CONCLUSION 

We have presented here a DFT study of the ground 
state structures of the undoped and Pd-doped rare-earth 
perovskite BaCeOs. The Wyckoff positions and space 
group assignment for BaCeOs are in agreement with ex- 
periment. The substitution of Ce by Pd without an 
accompanying O vacancy distorts the doped perovskite 
moreso than with an accompanying vacancy. When the 
vacancy is not adjacent to Pd, but is between two Ce, the 
structure is less stable. Located between a Pd-Ce pair, 
the vacancy causes the CeO^ square pyramid to tilt away 
from Pd, stabilizing a four-coordinate square-planar Pd 
geometry. The stabilization is about 0.95 eV relative to 
Structure 2. The coordination of Pd is also supported 
in the position of filled Pd (i-states just below the Fermi 
level of Structure 1 that are not filled in either Structures 
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FIG. 5: a) The HOMO of Structure 1, where the electron 
density is locahzed around the Pd cation, and its nearest Ce 
neighbors, b) The LUMO of Structure 1, composed solely of 
Ce /-orbit als. There is no electron density localized around 
Pd. c) An empty state above the LUMO, composed of Ce /- 
orbitals and a Pd d^2_y2 . In all figures, Pd is the green cation 
located at the right corner, Ce are white, Ba are light blue 
and O red. Bonds are omitted for clarity of the isosurface 
shown. 



2 or 3. This is also supported by the neutron pair distri- 
bution function data in Reference [17^ . The stabilization 
of oxygen vacancies by Pd doping leads to an enhanced 
rate of oxygen diffusion through the material [18 , which 
could be used in the design of a solid oxide fuel cell mem- 
brane or low temperature catalyst. 
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FIG. 6: Projected Density of States of the Pd d-states. Struc- 
ture 1 (top) shows that a filled band is being stabilized when 
an oxygen vacancy is placed between a Pd-Ce pair, as d-states 
above the Fermi level become less populated. This feature 
is not seen in either Structures 2 (middle) or 3 (bottom). 
In these two structures, there appear sharply intense bands 
around -5 eV, absent in Structure 1. These are filled states, 
lower in energy, that demonstrate the more ionic nature of 
the Pd in Structures 2 and 3. 
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